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1 Introduction 



We review efforts to unify both the Bardeen, Cooper & Schrieffer (BCS) and Bose-Einstein 
condensation (BEC) pictures of superconductivity. We have finally achieved this in terms of 
a "complete boson-fermion (BF) model" (CBFM) that reduces in special cases to all the main 
continuum (as opposed to "spin") statistical theories of superconductivity. Our BF model is 
"complete" in the sense that not only two-electron (2e) but also two-hole (2h) Cooper pairs 
(CPs) are allowed in arbitrary proportions. In contrast, BCS-Bogoliubov theory — which can 
also be considered as the theory of a mixture of kinematically independent electrons, 2e- and 
2h-CPs — allows only equal, 50%-50%, mixtures of the two kinds of CPs. This is obvious from 
the perfect symmetry about fi, the electron chemical potential, of the well-known Bogoliubov 
[1] v 2 (e) and u 2 (e) coefficients, where e is the electron energy. The CBFM is then applied to 
see: a) whether the BCS model interaction for the electron-phonon dynamical mechanism is 
sufficient to predict the unusually high values [2] of T c (in units of the Fermi temperature) of 
~ 0.01 — 0.1 exhibited by the so-called "exotic" superconductors [3] in both 2D and 3D — relative 
to the low values of < 10 -3 more or less correctly predicted by BCS theory for conventional, 
elemental superconductors; and b) whether it can at least suggest, if not explain, why "hole 
superconductors" have higher T c 's. 



Boson-fermion (BF) models of superconductivity as a BEC [4, 5] go back to the mid- 
1950's [6]-[9], pre-dating even the BCS-Bogoliubov theory [10, 11]. Although BCS theory only 
contemplates the presence of "Cooper correlations" of single-particle states, BF models [6]- [9], 
[12]-[20] posit the existence of actual bosonic CPs. Such pair charge carriers have been observed 
in magnetic flux quantization experiments with elemental [21, 22] as well as with cuprate [23] 
superconductors. But apparently no experiment has yet been done that distinguishes between 
electron and hole CPs. The fundamental drawback of early [6]- [9] BF models, which considered 
2e bosons in analogy with diatomic molecules in a classical gas mixture, is the notorious absence 
of an electron energy gap A. The gap first began to appear in later BF models [12]-[16]. With 
two [15, 16] exceptions, however, all BF models neglect the effect of hole CPs formulated on 
an equal footing with electron CPs to give us a complete BF model (CBFM) consisting of both 
bosonic CP species coexisting with unpaired electrons. 

Without going into a detailed justification we merely list several common "myths" in the 
theory of superconductivity that we tacitly disbelieve: 

1. With the electron-phonon dynamical mechanism transition temperatures T c < 45 K at 
most. For higher T c 's one needs magnons or excitons or plasmons or other electronic 
mechanisms. 

2. Cooper pairs (CPs): 

a) consist of negative-energy stable (i.e., stationary) bound states [24]; 

b) propagate in the Fermi sea with energy h 2 K 2 /2(2m) (Ref. [9], p. 94) where HK is 

the total or center-of-mass momentum (CMM) of the composite pair; 

c) have a linear dispersion E oc vphK, with vp the Fermi velocity, which is merely the 

acoustic mode in the ideal Fermi gas (+ interactions) with sound speed Vp/Vd in 
any dimensionality d; 

d) "...with K 7^ represent states with net current flow [25];" 

e) are not bosons (Ref. [26], p. 38). And most notoriously, that: 

3. Superconductivity is unrelated to BEC [27]. 

We question all of these assertions which will be discussed in greater detail elsewhere. 

2 The CBFM Hamiltonian 

The CBFM [15, 16] is described by H = Hq + Hi nt where the unperturbed Hamiltonian Hq 
corresponds to an ideal (i.e., noninteracting) gas mixture of fermions and both types of CPs, 
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two-electron (2e) and two-hole (2h), namely 

H = £ k 1 « 1 r 1 ,s 1 «k 1) s 1 +j2 E +( K ) b K b * - E e -W c k c k, (1) 

ki,si K K 

where K = ki + k 2 is the CMM wavevector, k = 1^ — k 2 ) being the relative one, while 
e k = h 2 k 2 /2m are the electron and E±(K) the 2e-/2h-CP energies. Here a£ s (aki,si) are 
creation (annihilation) operators for fermions and similarly 6 K (6k) and c K (ck) for 2e- and 
2h-CP bosons, respectively. 

Two-hole CPs are considered distinct and kinematically independent from 2e-CPs as their 
Bose commutation relations involve a relative sign change, in sharp contrast with electron or 
hole fermions whose Fermi anticommutation relations do not. In fact, holes have a dramatic 
effect even in the simple, elementary CP problem where they were originally neglected thereby 
giving [24] a negative-real-energj, stationary (i.e., infinite-lifetime) two-fermion bound-state. 
But if electrons and holes are treated equally and simultaneously through a Bethe-Salpeter 
(BS) equation (see, e.g., Ref. [28] p. 131) in the ideal Fermi gas (IFG) ground-state about 
which the CPs are defined, the resulting energy is pure imaginary [11], [29] — implying an obvious 
instability. The IFG-based CP problem is thus meaningless if particles are taken on an equal 
footing with holes, as consistency would demand. However, a similar BS treatment not about 
the IFG but about the BCS ground-state yields [30] real (but positive, as with a "quasi-bound 
state in the continuum") 2e- and 2h-CP energies, along with an imaginary part that is nonzero 
only for K ^ signifying a finite lifetime, but zero for K = implying permanent pairs. Thus, 
the CP problem is vindicated in a very natural, physical way via the BS equation. 

The interaction Hamiltonian ifj nt consists of four distinct interaction vertices, each with 
two-fermion/one-boson creation or annihilation operators, depicting how unpaired electrons 
(subindex +) [or holes (subindex — )] combine to form the 2e- (and 2h-CPs) assumed in the 
system of size L, namely 

H int = L~ V2 /+( A; ){ a k+lK,T a -k+|K,| 6K + a -k+iKa a k+|K, T & K} 
k,K 

+L- 3 / 2 J2 /-WK+iK^-k+iK^K + a -k+|K,l a k+|K,T C K}- (2) 
k,K 

Note that the fermion-pair interaction H int is reminiscent of the Frohlich (or Dirac QED) 
interaction Hamiltonian involving two fermion and one boson operators but with two types of 
CPs instead of phonons (or photons). But in contrast with Frohlich and Dirac there is no a 
conservation law for the number of unpaired electrons, i.e., [H int , Yl ^ki^ s Oki.sJ 7^ 0. (Note 

ki,si 

too that [H int , kia£ a kl , Sl ] = and [H int , Yl s i a kl , Sl a ki,sJ = 0.) Just as the Frohlich 

ki,si ki,si 

(or Dirac) interaction Hamiltonians are the most natural ones to use in a many-electron/phonon 
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(or photon) system, one can conjecture the same of (2) for the BF system under study. Indeed, 
this Hi n t has formally already been employed under various guises by several authors [12], [31]- 
[34]. The energy form factors f±(k) are essentially the Fourier transforms of the 2e- and 2h-CP 
intrinsic wavef unctions, respectively, in the relative coordinate between the paired fermions of 
the CP. Here they are taken simply as 

f+{e) = jf torE f <e<E f + 5e, 
' otherwise, 

f-(e) = { f for ^/- & <^<^ (4) 
' otherwise, 

with Ef and 5e phenomenological dynamical energy parameters (in addition to the positive 
coupling parameter /) related to the bosonic CPs through Ef = ^[E + (0) + £L(0)] and Se = 
7}[E + (0) — £L(0)], where E±(0) are the (empirically unknown) zero-CMM energies of the 2e- 
and 2h-CPs, respectively. Clearly £±(0) = 2Ef±5e. The quantity Ef will serve as a convenient 
energy scale and is not to be confused with the Fermi energy Ep — \mv F = IvbT f where T F 
is the Fermi temperature. The Fermi energy E F equals 7rh 2 n/m in 2D and (ft 2 /2m)(37r 2 n) 2//3 
in 3D, with n the total number-density of charge- carrier electrons. The quantities Ef and E F 
coincide only when perfect 2e/2h-CP symmetry holds. 



3 Main statistical theories as special cases of CBFM 

The interaction Hamiltonian (2) can be further simplified by keeping only the K = terms. One 
then applies the Bogoliubov "recipe" (see, e.g., [28] p. 199) of replacing in the full hamiltonian 
H = H + H int all zero CMM creation and annihilation operators by c-numbers: y/No and 
v/Mq for 2e- and 2h-CP operators, where N (T) and M (T) are the number of zero-CMM 2e- 
and 2h-CPs, respectively. Minimizing with respect to the independent variables N and M the 
so-called thermodynamic (or grand) potential associated with the full Hamiltonian H , as well 
as keeping the total number of electrons fixed and thereby introducing the electron chemical 
potential /z, yields a set of three coupled, transcendental, integral equations (Ref. [15], Eqs. 
7-8). These three equations embody the CBFM. Two of these are coupled gap-like equations 
involving the 2e-CP and 2h-CP BE-condensed boson number densities n (T) = N (T)/L d and 
mo(T) = M (T)/L d , linked together through an electron energy gap A. The third equation can 
be cast as a number equation of the form 2ne(T) — 2m#(T) +ny(T) = n involving both 2e and 
2h boson number-densities but now for all energy states, where nf(T) is the number-density of 
unpaired electrons. Most significantly, these three equations contain five different theories as 
special cases, see flow chart in Fig. 1. For perfect 2e/2h CP symmetry ub(T) = itlb{T) implies 
[15] that no(T) = m (T) and that Ef coincides with /i. The CBFM then reduces to: i) the gap 
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and number equations of the BCS-Bose crossover picture [35] for the BCS model interaction — if 
the BCS parameters V and Debye energy Hud are properly identified with the CBFM dynamical 
parameters through f 2 /25e and 5e, respectively The crossover picture for unknowns A(T) and 
n(T) is now supplemented by the key relation A(T C ) = fy^n (T c ) = f \Jm^(T\). If in addition 
one imposes that /i ~ E F , as occurs for weak coupling from the number equation, the crossover 
picture is well-known to reduce to: ii) ordinary BCS theory. Thus, the BCS condensate is 
precisely a BE condensate whenever ne(T c ) = rriB(T c ) and the coupling is weak. 

On the other hand, for no 2h-CPs present the CBFM reduces [15] also to: iii) the BEC 
BF model in 3D of Friedberg and Lee [13, 14] characterized by the relation A(T) = f^n (T). 
With just one adjustable parameter (the ratio of perpendicular to planar boson masses) this 
theory fitted [14] cuprate T C /T F data quite well. When / = the CBFM reduces to both: iv) 
the ideal BF model of Ref. [19, 20] that predicts nonzero BEC T c 's even in 2D, as well as to: v) 
the familiar T c -formula of ordinary BEC in 3D, albeit as an implicit equation with the boson 
number-density being T-dependent. Figure 2 displays the T c prediction [20] in 2D for cuprate 
superconductors of special case (iv), with no adjustable parameters. 



4 BEC limit of all electrons paired 

The general BEC T c -formula for noninteracting bosons in ^-dimensions of energy e K = C S K S , 
s > 0, is [18] 



C s 

T r = T ± 

k 



- c 



B 



sT(d/2)(2ir) d 
2^T{d/s)g d/s {l) nB , 



s/d 

(5) 



where % is the boson number- density and the Bose integral [36] g a (z) = ^ z l /l a , with z = 

i=i 

e p.B/k B T j g foe "fugacity" and /i# the boson chemical potential. For z = 1, <?<r(l) = the 
Riemann Zeta-function, if a > 1, while for < o < 1 the infinite series g a {l) diverges. Eq. (5) 
is formally valid for all d > and s > 0. Hence, for < d < s, T c = since ^/ s (l) = oo for 
d/s < 1 but T c is otherwise finite. We stress that as a consequence of the former all 2D T c 
predictions in Fig. 2 (except the BCS one that survives for all d > 0) would collapse to zero 
had s = 2 been used in 2D instead of the correct s — 1. For s = 2, Ci = h 2 /2mB (5) leads 
to the familiar 3D result T c ~ 3.31h 2 n 2 ^ 3 jm^kB since £(3/2) ~ 2.612. Recalling that ksTp = 
h 2 k 2 F /2m with k F = [2 d ~ 2 n d ^ 2 dT(d/2)n] 1 ^ d , then if m B = 2m and ris = n/2 (all electrons 
paired) for s = 2 (5) gives T C /T F = \\2j 1 dT(d)g d/2 (l)] 2 / d = for d < 2 since g d/2 (l) = oo for 
d/2 < 1. For d = 3 we arrive at another familiar result T C /T F = \ [2/3r(3/2)C(3/2)] ~ 0.218 
(see dashed line in "Uemura plot" of Ref. [2], Fig. 2). This value appears marked in Fig. 4 as 
a black triangle. 

We now focus on s = 1. For the boson energy rj to be used below the leading term in the 
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many-body Bethe-Salpeter (BS) CP dispersion relation is linear, i.e., r\ ~ (X/2n)hv F K [see Ref. 
[30] for the derivation in 3D which gives r] ~ (\/A)hv F K\. Here A = VN(E F ) where N(E F ) is 
the electron density of states (DOS) (for one spin) at the Fermi surface. Note that the boson 
energy r\ is linear in CMM K — and not the quadratic h 2 K 2 /4m appropriate for a composite 
boson of mass 2m moving not in the Fermi sea but in vacuum [6]-[9], [13]-[16]. The quadratic 
holds only when Ep is strictly zero [17], i.e., when no Fermi sea is present. These linearly- 
dispersive CPs are commonly confused with the also linearly-dispersive sound phonons of the 
collective excitation sometimes referred to as the Anderson-Bogoliubov-Higgs (ABH) (Ref. [11] 
Sec. 3; [37, 38]) mode, which for zero coupling reduces [39] to the IFG result hv F K/\fd. The 
IFG sound speed c = v F /\fd also follows directly from the zero-temperature IFG pressure 
P = n 2 [d(E /N)/dn] = 2nE F /(d + 2) via the familiar thermodynamic relation dP/dn = mc 2 , 
where E is the ground-state energy and as before n = N/L d = k d F /2 d - 2 n d / 2 d T(d/2) is the 
fermion- number density. But the above results with rj oc Xhv F K in fact refer to actual "moving" 
(or "excited") CPs in the Fermi sea. Both kinds of distinct soundwave-like solutions — moving 
CPs and ABH phonons — appear in the many-body BS ladder-summation scheme of Ref. [30]. 
Note also that the BS CP linear dispersion coefficient X/2tt in 2D (or A/4 in 3D) contrasts 
markedly with the coupling-independent 2/tt coefficient in 2D (or 1/2 in 3D, as first quoted in 
Ref. [26], p. 33) obtained [40] in the simple CP problem [24] which ignores holes. Thus (5) 
again with ub = n/2, for s = 1 and C\ = Xb(d)hv F with 6(2) = l/2n and 6(3) = 1/4, yields 
T C /T F = 2\b(d)/[dT(d)((d)] 1/d which (for A = 1/2) is ~ 0.088 if d = 2 since ((2) = tt 2 /6, and 
~ 0.129 if d = 3 since £(3) ~ 1.202. These two values for T C /T F will appear as the uppermost 
black squares in Fig. 4 marking the BEC limiting values if all electrons in our 2D or 3D many- 
electron system were imagined paired into noninteracting bosons formed with the BCS model 
interelectron interaction. For A = 1/4 the lowermost black squares in Fig. 4 apply. 

5 Enhanced T c 's from the CBFM 

We now apply this very general CBFM to exhibit the sizeable enhancements in T c s over BCS 
theory that emerge for moderate departures from perfect 2e/2h-pair symmetry for the same 
interaction model. The pair-fermion interaction (2) with (3) and (4) bears a one-to-one corre- 
spondence with the more familiar "direct" interfermion electron-phonon interaction, mimicked, 
e.g., in the BCS model interaction (whose double Fourier transform is a negative constant — V 
nonzero only within an energy shell 2Huj Fj about the Fermi surface, with up, the Debye fre- 
quency) if [15, 16] we set f 2 /28e = V and Se = fkoo- The familiar dimensionless BCS model 
interaction parameters A = N(E F )V and hwrj/E F are then recovered. 

The three coupled equations of the CBFM determining the <i-dimensional BE-condensate 
number-densities n (T) and m (T) of 2e- and 2h-CPs, respectively, as well as the fermion 
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chemical potential «(T), were solved numerically in 3D for A = 1/5 and ftu B / ' E F = 0.001 in 
Ref. [16] assuming a quadratic boson dispersion relation rj = ft 2 K 2 /4m. For this case Figure 3 
maps the phase diagram in the vicinity of the BCS T c value (marked BCS-B in the figure) at 
An = n/rif — 1 = (corresponding to perfect 2e/2h-CP symmetry) where rif is the number of 
unpaired electrons for zero gap and zero temperature, provided that rif < n [15, 16]. Besides the 
normal phase (n) consisting of the ideal BF gas described by H , three different types of stable 
(plus several metastable, i.e., of higher Helmholtz free energy) BEC-like phases emerged. These 
are two pure phases of either 2e- (s+) or 2h-CP (s— ) BE-condensates, and a lower temperature 
mixed phase (ss) with arbitrary proportions of 2e- and 2h-CPs. Of greater physical interest are 
the two higher-T c pure phases so that we focus below only on them. For each pure phase at a 
critical temperature we have either A(T CS+ ) = / ^/no(T cs+ ) = or A(T CS _) = /^/mo(T cs _) = 0, 
where A(T) is the electronic (BCS-like) energy gap. Their intersection gives the BCS T c value 
of 7.64 x 10'% where k B T f = E f = (ft 2 /2m)(3n 2 n f ) 2 / 3 . 



5.1 Two dimensions (2D) 

In 2D the one-spin electronic density of states (DOS) is constant, namely N(e) = mjl-nft} . Us- 
ing the aforementioned BS CP linear dispersion relation rj ~ (\/27t)HvfK we get for the bosonic 
DOS M(rj) = (l/2ir)K(dK/dr)) ~ (2n/ \ 2 ft 2 v 2 F )r] instead of the constant that follows in 2D from 
quadratic dispersion. Employing Ef = nft 2 nf/m = k B Tf as energy/density/temperature scal- 
ing factors, and the relation n/rif = (Ep/ Ef) d ^ 2 to convert quantities such as T c /Tf to T c /Tp, 
where E F = k B T F , the two working equations for the pure 2e-CP phase [i.e., m B (T c ) = 0] be- 
come (with all quantities dimensionless, energies in units of Ef and electron particle-densities 
in units of nf) 

1 + ftuj D /2 - u = \(ftw D /2) I dx- — - — -tanh^— - (6) 

J \x - u| 2T C 
i 

oo oo 

1 1 Ml- tanh^] + / fa [coth E + jOj^gZiziO _ !] = ,, (7) 



These are just the gap-like equation associated with 2e-CPs and its corresponding number 
equation. 

For the pure 2h-CP phase (i.e., n B = 0) the two working equations are 



i 



u - 1 + ftw D /2 = X(ftuj D /2) f dx- — - — -tanh ^ ^ , (8) 

J \X fi\ 21 c 



1 — Hu>£, 
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oo oo 

1 f , r „ . x — u, 7r 2 1 /" , r , x — 2(1 — hjJr>/2 — a) . . 

- y cfa[l - tanh^] - — - y ^[coth 1 2Tc D/ ^ - 1] = n. (9) 



2 

o o 



which are the gap-like equation associated with 2h-CPs and its corresponding number equation. 
Note that (7) and (9) are quadratic in n, and that all integrals there are exact, namely 



oo 

/ 



dx[l - tanh^-— ^] = /i + 2T c ln[2cosh( / u/2T c )], (10) 
z i c 



and 



J dxx[coth^- - 1] = 2T 2 c g 2 {e~^) (11) 
o 

where the Bose function g a (z) was defined after (5); it is designated in Ref. [41] as PolyLog[a, z\. 
The integrals in (6) and (8) were performed numerically. In 2D we use the two extreme values 
of A = 1/4 (lower set of curves in Fig. 4) and = 1/2 (upper set of curves), and tko D /E F = 0.05 
(a typical value for cuprates), to compute from equations (6) to (9) the T c /Tp vs. n/rif phase 
diagram which is graphed in the figure (left panel for 2D) for both 2e-CP (dashed curve) 
and 2h-CP (full curve) pure, stable BEC-like phases. The value n/rif = 1 corresponds to 
perfect 2e/2h-CP symmetry. The T c value where both curves n (T c ) = m (T c ) = intersect is 
marked by the large dots in the figure; these values are consistent with those gotten from the 
familiar BCS expression T C /T F ~ 1.13A(hv D / E F ) exp(-l/A) ~ 0.001 for A = 1/4, and 0.008 for 
A = 1/2, for hwo/Ep = 0.05. [Thes values differ little from those from the exact BCS (implicit) 
T c -formula (Ref. [28], p. 447) 1 = A J^ D / 2kBn dxx' 1 tanhx.] Cuprate data empirically [42] 
fall within the range T C /T F ~ 0.03 — 0.09. Thus, moderate departures from perfect 2e/2h-CP 
symmetry enable the CBFM to reach quasi-2D cuprate empirical T c values, and quite likely also 
room temperature superconductivity, without abandoning electron-phonon dynamics — contrary 
to popular belief. Compelling evidence for a strong, if not sole, phonon dynamical component 
in cuprates has recently been reported [43] from angle-resolved-photoemission data. 



5.2 Three dimensions (3D) 

In3DiV(£) = (l/27r 2 )k 2 (dk/de) = (mi /2^-ir 2 h 3 ) y/e and, analogously as before, Ef = (h 2 /2m)x 
(37r 2 n/) 2 / 3 = ksTf which again differs from Ep = (h 2 /2m)(3'K 2 n) 2 ^ 3 = ksTp except when per- 
fect 2e/2h-CP symmetry holds when they coincide, whereas the leading term in the BS CP 
boson dispersion energy is now the linear expression rj ~ (\/4)HvfK [30] so that M(i]) = 
(l/2TT 2 )K 2 (dK/dr]) ~ (32/7r 2 X 3 h 3 v F )r] 2 . The above working equations for the pure 2e-CP 
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phase now in 3D (all quantities again dimensionless) are 

If 1 II 

l + hoj D /2- fi = \{hjj D /2)— - / dxy/x-. rtanh ^"^ 1 , (12) 

n ' J \x — fi\ 2T C 

1 

oo oo 

- y rfxv^[l - tanh^] + -- y ^ 2 [coth 1 ^ ^ - 1] = n, (13) 



while for the pnre 2h-CP phase they are 

i 

1 f X 1 3? I 

H-l + hu) D /2 = \{hw D /2)^r K I dxy/x-. -tanh F "l^ 1 , (14) 

n 1 ' 6 J \x — ji\ 2T C 



oo oo 



3 f , /— r , x — a 1 121/", 9r , x — 2(1 — Hojd/2 — u) 

- y dxv^[l - tanh^] - -- y ^ 2 [coth 1 ^ ^ - 1] = n. (15) 



Results in 3D are reported only for the 2e-CP BEC case for the same extreme values of A = 1/4 
and 1/2 as in 2D but now for Hujd/E f = 0.005. In Fig. 4 for 3D (right panel) the dashed 
curves are the 2e-CP BEC phase boundaries. The large dot again marks the BCS T C /T F values 
of 0.0001 for A = 1/4 and 0.0008 for A = 1/2. Empirical data for both exotic and conventional, 
elemental superconductors in 3D are taken from Ref . [2] . We see that whereas BCS theory can 
roughly reproduce T c jTp values well for the latter, it takes moderate departures from perfect 
2e/2h-CP symmetry to access 3D exotic superconductor T c /Tp values, which empirically [2] 
fall within the range ~ 0.01 — 0.1. This is much larger than the range < 0.001 for conventional 
(elemental) superconductors, also shaded in the right panel of the figure. 



6 Hole superconductivity 

Finally, we address the unique but mysterious role played by holes in superconductors in gen- 
eral. For example: a) of the cuprates those that are hole-doped have transition temperatures 
T c about six times higher than electron-doped ones; and b) in fullerite (an fee crystal of C 60 
fullerenes) T c is now claimed to be more than three times higher with hole rather than electron 
doping, as recently observed [44] with the so-called "field-effect transistor" technique of inject- 
ing holes. And even in conventional superconductors [45] c) over 80% of all superconducting 
elements have positive Hall coefficients (meaning hole charge carriers); and d) over 90% of non- 
superconducting metallic, non-magnetic elements have electron charge carriers. This greater 
"efficiency" of holes in producing higher T c 's is clearly predicted in Fig. 4 at least for 2D 
superconductors, at least insofar as 2h-CP BE condensates exhibit higher T c 's than 2e-CP 
ones. 
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7 Conclusions 



This brief review has sketched how five statistical continuum theories of superconductivity- 
including both the BCS and BEC theories — are contained as special cases in a single theory, the 
"complete boson-fermion model" (CBFM). This model includes, for the first time, both two- 
electron and two-hole pairs in freely variable proportions, along with unpaired electrons, all in 
chemical/thermal equilibrium. The BCS condensate (characterized by a single equation, namely 
the T-dependent gap equation) follows directly as a BE condensate through the condition for 
phase equilibria when both 2e and 2h pair numbers are equal at a given temperature and 
coupling — provided the coupling is weak enough such that the electron chemical potential is 
roughly the Fermi energy. Ordinary BEC theory, on the other hand, is recovered from the 
CBFM when hole pairs are neglected, fermion-pair coupling is made to vanish, and the limit 
of all electrons being paired into bosons is taken. 

The practical outcome of this BCS-BEC unification via the CBFM is then threefold: a) 
enhancements in T c , by more than an order-of-magnitude in 2D, and more than two orders- 
of-magnitude in 3D, are obtained for the same electron-phonon dynamics mimicked by the 
BCS model interaction — provided only that one departs moderately from the perfect 2e/2h- 
pair symmetry to which BCS theory is intrinsically restricted; b) these enhancements in T c 
fall within empirical ranges for 2D and 3D "exotic" superconductors, whereas BCS T c values 
continue to lie low and within the empirical ranges for conventional, elemental superconductors; 
and c) hole-doped superconductors are predicted to have higher T c 's than electron- doped ones, 
in agreement with observation. 
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Figure 1: Flow chart of how the CBFM reduces in special cases to the statistical, continuum 
models of superconductivity discussed in text, thereby displaying how both BCS and BEC 
theories can be unified. 
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Figure 2: Critical 2D BEC-like temperature T c in units of Tp for the BCS model interaction 
with A = 1/2 for varying hwn/Ep = Op/Tp for: the pure unbreakable-boson gas with some 
and with all fermions paired; for the breakable-boson gas; and for the boson-fermion mixture 
(thick full curve labeled "binary gas") in thermal/chemical equilibrium, all as described in Ref. 
[20] for the original (simple) CPs where C\ = (2/ir)hvF- Dashed curve is the BCS theory T c , 
and cuprate experimental data are taken from Ref. [42]. 
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Figure 3: Phase diagram [16] with 3D superconducting critical temperatures T cs+ , T cs _, T css+ , 
and T css _ as functions of An = n/rif — 1 as defined in text, in the vicinity of the BCS T c value, 
assuming a quadratic boson dispersion, for A = 1/5 and Hud/Ef = 0.001. 
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Figure 4: Phase diagrams in 2D and 3D for temperature (in units of T F ) and electron density 
(in units of nj as defined in text) showing the phase boundaries of T c 's for the pure 2e-CP 
BEC phases (dashed curves) determined by A(T C ) = f \/n^{T^j = 0, and the pure 2h-CP BEC 
phase (in 2D only) given by A(T C ) = fy/m (T c ) = for A = 1/4 and 1/2 with hw D /E F = 0.05 
in 2D and 0.005 in 3D. Intersections corresponding to no(T c ) = mo(T c ) giving the BCS T c 
approximately are marked by black dots, while black squares mark the BEC limit where all 
electrons are imagined paired into 2e-CP bosons, and the black triangle marks the familiar 3D 
limit as determined in Sec. 4. 
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